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ABSTRACT 

We present the first results of the Herschel open time key program, Galactic Observations of Terahertz C* (GOT C+) survey of the 
[CIl] 'P3/2-^Pi/2 fine-structure line at 1.9 THz (158 yum) using the HlFl instrument on Herschel. We detected 146 interstellar clouds 
along sixteen lines-of-sight towards the inner Galaxy. We also acquired HI and CO isotopologue data along each line-of~sight for 
analysis of the physical conditions in these clouds. Here we analyze 29 diffuse clouds (Ay < 1.3 mag. ) in this sample characterized by 
having [ClI] and HI emission, but no detectable CO. We find that [CIl] emission is generally stronger than expected for diffuse atomic 
clouds, and in a number of sources is much stronger than anticipated based on their HI column density. We show that excess [CIl] 
emission in these clouds is best explained by the presence of a significant diffuse warm Ht, dark gas, component. This first [CIl] 158 
jum detection of warm dark gas demonstrates the value of this tracer for mapping this gas throughout the Milky Way and in galaxies. 
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1. Introduction 

Interstellar gas plays a crucial role in the life cycle of Galaxies 
providing the material for star formation and as a repository of 
gas ejected by stars as part of their evolution. The diffuse atomic 
gas has been mapped in HI 21 cm surveys, and dense molecu- 
lar H2 clouds have been mapped indirectly with molecular trac- 
ers, primarily '-CO (c.f. Dame et al. 2001) and in more limited 
surveys in its isotopologues. Missing, however, is a widespread 
spectral tracer of H2 not located in regions with conditions ap- 
propriate to form CO, the so-called "dark gas" (Grenier et al. 
2005). One of the best candidates for this tracer is C^, which is 
widely distributed in the Galaxy as shown by COBE FIRAS all 
sky (Bennett et al. 1994) and BICE inner Galaxy (Nakagawa et 
al. 1998) observations of its 158 //m ^P3/2-^Pi/2 fine structure 
line. In addition, [CIl] is a density- and temperature-sensitive 
probe of diffuse clouds and photon dominated regions (PDRs). 
While COBE revealed that this line is the brightest far-IR line 
in the Galaxy, low spectral (Av w 1000 km s"') and spatial (7°) 
resolution could not reveal individual cloud components. BICE 
data, while slightly better in this respect (175 km s"' and 15'), 
is not adequate to resolve clouds along the line-of-sight (LOS). 
Heterodyne receivers are critical to provide the necessary ve- 
locity resolution to identify and study individual clouds along 
the LOS; however, prior to operation of the HIFI instrument (de 
Graauw et al. 2010) on the Herschel (Pilbratt et al. 2010) only a 
handful of high spectral resolution [CIl] spectra were available 
from bright HII regions (c.f. Boreiko & Betz 1991). 

Here we report the first results of a large-scale survey of 
[CIl] 1.9 THz line emission in the Galaxy being conducted un- 
der the Herschel open time key program. Galactic Observations 
of Terahertz C+ (GOT C+). The full GOT C+ program will 
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observe about 900 LOS in the Galaxy, divided into four sub- 
programs (Langer et al(2010): 1) a Galactic disk uniform vol- 
ume sampling in longitude (covering all 360°), at b = 0, +0.5°, 
±1°, and ±2°; 2) strip maps in the inner Galaxy; 3) a sampling 
of Heiles high-latitude clouds studied in absorption; and, 4) the 
Galactic warp in the outer Galaxy. To date we have obtained high 
spectral resolution spectra of [CIl] towards 16 LOS towards the 
inner Galaxy near longitude 340° (5 LOS) and 20° (11 LOS) 
taken during the performance verification phase and priority sci- 
ence program. In this small sample (< 2% of the complete GOT 
C-l- survey) we detected 146 clouds in [CIl] emission. These 
data cover four broad categories of clouds as defined by their HI, 
[CIl], and CO signatures: a) diffuse atomic clouds with high HI 
to [CIl] intensity ratio and no detectable CO; b) diffuse molecu- 
lar clouds with relatively low HI to [CIl] intensity ratio (due to 
relatively strong [CIl]) and no '^CO emission (Av< 1.5 mag. ); 
c) transition molecular clouds and PDRs detected in HI, [CIl], 
and '^CO, but not '^CO; and, d) dense molecular clouds detected 
with HI, [CIl], '^CO, '^CO, and sometimes C'**0. 

In this Letter we analyze 29 diffuse clouds (Av<1.3 mag.) 
characterized by the presence of [CIl] and HI emission, but hav- 
ing no CO emission. We find that [CIl] emission is generally 
stronger than expected for diffuse atomic clouds (T ~ 100 - 
200 K, n ~ 100- 200 cm"^), and in about one-third of the sources 
[CIl] emission is much stronger than anticipated. This [CIl] 
emission indicates the presence of warm "dark gas" (H2 with- 
out CO) associated with HI. Velusamy et al. (2010) discuss the 
transition C^-'^CO clouds and their "dark gas" envelopes, and 
Pineda et al. (2010) discuss the dense molecular cloud PDRs. 



2. Observations 

The 16 LOS observations of [CIl] at 1.9 THz were made 
with Herschel HIFI band 7b using Load Chop (HPOINT) with 
SkyRef offset by 2° off the plane, for better cancelation of instru- 
mental spectral baselines. The data were reduced with HIPE ver- 
sion 3 and a fringe-fitting tool to remove standing waves, and the 
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Fig. 1. [CII] spectra obtained with Herschel HIFI for the GOT C+ program at /=345.65° and b=0°, along with the CO data from 
the Mopra telescope and HI surveys. Many different types of interstellar clouds can be seen in this LOS; diffuse clouds with HI and 
C^ emission, but without, or barely any, CO; strong [CII] emission towards dense molecular clouds, as indicated by the presence of 
'^CO and in some cases C"*0. The Gaussian decomposition for [CII] is also shown in the upper panel. 



h- and v-polarizations combined where available. All [CII] data 
are corrected for the antenna's main beam efficiency of 0.63. We 
used the wide band spectrometer (WBS) with 0.22 km s~' at 1.9 
THz. The 16 LOS are: 1) b=0 at longitudes: 337.82°, 343.04°, 
343.91°, 344.78°, 345.65°, 18.26°, 22.60°, 23.47°, 24.34°; and, 
2) b = 0.5° at 24.34°, b=-0.5° at 18.26°, 23.47°, b=+1.0° at 
22.60°, 24.34°, and b=-1.0° at 18.26°, 23.47°. Further details 
are in Velusamy et al. (2010). We also observed the 7 = 1 — > 
transitions of '^CO, '^^CO, and C'^O toward each LOS with the 
ATNF Mopra 22-m Telescope (details in Pineda et al. 2010), 
with an angular resolution of 33". We obtained HI data from 
public sources (McClure-Griffiths et al. 2005, Stil et al. 2006). 

A typical example of our data is shown in Figure 1 for 
1=345.6522°, b=0.0°, where we plot the [CII] antenna tempera- 
ture (corrected for main beam efficiency) along with those of HI, 
and the CO isotopes (other examples are shown in Pineda et al. 
2010 and Velusamy et al. 2010). This LOS passes through the 
edge of the bar in the inner Galaxy. The [CII] features are clearly 
blends of many cloud components. It can be seen in Figure 1 that 
a wide variety of clouds are detected in [CII], for example: be- 
tween V/jr - -20 and 10 km s"' there is strong [CII] emission 
along with CO and '^CO, whereas between -100 and -140 km 
s"' four narrow [CII] features appear which have weak '^CO as- 
sociated with them, while between -50 and -90 km s"' there are 
five features with only [CII] and HI, and no detectable '^CO. 

3. Results and analysis 

We identified our clouds from Gaussian decompositions (see 
Figure 1) of [CII] and '^CO (where available) and use these 
to characterize their velocities, Vlsr^ and linewidth Av (see 
Velusamy et al. 2010 for details). With this approach we de- 
tected [CII] in 146 components at the 3-cr level or better We 
did not fit the HI spectral profiles, but obtained the HI intensi- 
ties in each cloud for comparison with other lines, by integrating 
within the velocity width (Av) centered at their Vlsr defined by 
the [CII] (and ^^CO) fines. For each feature we derive the line pa- 
rameters, Tpeaki^, Av(km s"'), I = /Tdv (units of K km s"')), 
as well as those for HI and CO (where detected). Here we focus 
on 29 out of 35 [CII] components that do not have any known 
molecular gas as traced by '^CO (the other six have marginal 



[CII] intensities, or other issues). The linewidths (FWHM) for 
these 29 diffuse C^ clouds, have a mean value of 3.4 km s"' , and 
range from 1.4 to 5.5 km s~^ In Figure 2 we plot the integrated 
intensity I(CII) versus I(HI) for all 29 diffuse cloud features. It 
can be seen that [CII] is not strongly correlated with I(HI). This 
result is not what one would expect if all the C^ is only in HI 
clouds. Furthermore, there is a large degree of scatter, and, sev- 
eral sources have large [CII] emission at small N(HI), which we 
suggest is best explained by "dark gas" - H2 without CO. 

The identification of "dark gas" in our [CII] data can be un- 
derstood by using simple models of the radiative transfer to re- 
late the HI and [CII] intensities, and then consider the effects 
of adding an additional H2 layer containing C^. We calculate 
the column density of ionized carbon assuming the gas is only 
atomic, from the relation, N(C+) = X(C+) N(HI), where X(C+) 
is the fractional abundance of C^ in the gas with respect to 
HI. Observations in the local Interstellar Medium (ISM) lead to 
X(C+) in the range 1.4 to 1.8 x lO^^^ (Sofia et al. 1997), and 
here we adopt a value of X(C^)=1.5 x 10 **. The column den- 
sity of HI can be estimated in the optically thin limit, where the 
brightness temperature is less than the kinetic temperature, as 
A^(HI) = 1.82xl0i^/(HI)(cm-2). We can also directly calculate 
the N(C^) column density from I(CII) if we know the temper- 
ature and density of the gas. If the gas arises only from the HI 
cloud the N(C^) solutions should be consistent within a reason- 
able range of density and temperature. We derived the column 
density for N(C^) in the optically thin limit (the effects of opac- 
ity are discussed below) with no background radiation field (the 
CMB is negligible). 



A^„(C+) = 



hc^Au, 



/ 



T„t(CII)dv (cm"2) 



(1) 



where T„,fc is the main beam antenna temperature, N„ is the upper 
level column density. The total N(C^) is derived by solving the 
excitation of the ^¥3/2 state (in the optically thin limit). 



N(C^) = 2.9 X lO'^Xl + 0.5(1 + ncr/n)e 



i&E/kT) 



)I(CII)(cm-^) (2) 



The critical density (for HI or H2) is defined as, Wci = 
^i(// (o'v)„/cm"-', where A„/ = 2.4x10"^ s"' and cr is the col- 
lisional de-excitation cross-section for HI or H2, respectively. 
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For HI collisions with C^ we used the (crv)„; calculations of 
Barinovs et al. (2005) and derived a fit, 3.8x10^1" T°i^ cm^ s"', 
good to 20% over the temperature range 50 - 300K, yielding 
n„{HI) = 3.1 X 10^ cm^^ at lOOK. The critical density forH2 + 
C+ de-excitation is twice that of HI (Flower 1988, 1990). 

In Figure 2 we plot the predicted I(CII) versus I(HI) as a 
function of temperature and density, assuming that N(C^) is 
given by, N(C+) = X(C+)N(HI). The solid Unes show results for 
typical diffuse atomic cloud conditions, T = 100 to 300K and n 
= 100 to 300 cm^l The lowest curve is for (T,n) = (lOOK, 100 
cm"^), a likely lower bound for these diffuse atomic hydrogen 
clouds; all of our observed I(CII) fall above this line. The result 
for somewhat denser diffuse clouds, (T,n) = (lOOK, 300 cm"^), 
is shown by the middle solid line in Figure 2; about 25% of our 
sources fall along this line. We also plot results for a very warm 
HI cloud at (T,n) - (300K, 300 cm"^) (upper solid line) to give 
an idea of what physical conditions might be needed to explain 
a larger percentage of our [CII] clouds. About half the observed 
[CII] clouds fall below this line. However, this solution requires 
a rather hot diffuse atomic cloud, with a pressure, nT = 9x10"* 
that is much greater than the generally accepted value. At best it 
explains about half the [CII] sources as arising from atomic gas, 
but it is probably an unrealistically high combination of (n,T); 
and, the remaining sources require [CII] from additional cloud 
components. In summary, [CII] emission from diffuse HI gas 
alone does not explain that observed in much of our sample. 

The analysis above assumes optically thin HI and [CII] emis- 
sion. We estimate the opacity correction for the [CII] lines by 
starting with the optically thin derived N(C^) and using the 
RADEX code (van der Tak et al. 2007) to calculate t(CII) for a 
range of diffuse cloud densities and temperatures. For I(CII) < 5 
K km s" ' , T < 0.5 under most density and temperature conditions 
for a typical linewidth of 3 to 4 km s"', and the escape proba- 
bility for [CII] photons is > 80%. Thus, any corrections to the 
column density are < 20%. The HI opacity is t(HI) =5.5x10"''' 
N(HI)/(T^ Av), or, using the optically thin estimate for N(HI) 
above, t(HI) =I(HI)/(TxAv), where T^ is the excitation temper- 
ature. We do not know excitation temperature for HI, but it is 
reasonable to set it to 100 - 150K, which yields t(HI) < 0.5, 
for I(HI) < 300 K km s"' , and a maximum of 0.8 for the sources 
near I(HI) = 500 Kkms '. In summary, these corrections are 
small, and maintain approximately the same N(C^)/N(HI) ratio 
as the optically thin case, which does not resolve the issue of ex- 
cess [CII] emission within the context of purely atomic diffuse 
clouds. 

Unless we adopt a very high density, n > 10^ cm"^ and/or 
high temperatures (T > 300K), some other gas contribution is 
necessary to explain the [CII] emission. The clouds with I(CII) 
> 2.5 Kkms ' and I(HII) < 300 Kkms ' cannot even be ex- 
plained by invoking such high HI densities and temperatures. 
Furthermore, we would be faced with explaining how such dense 
clouds, not associated with HII regions, could have such high 
temperatures. We suggest that [CII] in these stronger C^ emit- 
ting clouds, is tracing a significant component of warm molec- 
ular hydrogen gas, or warm "dark gas". Early cloud-chemical 
models predicted such an H2 gas layer with C^ and little or no 
CO (see Langer 1977 and references therein), although the term 
"dark gas" was not invoked until much later. Since then, ever 
more accurate and sophisticated cloud models have been devel- 
oped, but the basic structure of the HI to H2 and C^ to CO-CI 
profile remains the same (e.g. Tielens & Hollenbach 1985, van 
Dishoeck & Black 1988, Leeet al. 1996, Wolfire et al. 2010). 

To estimate the amount of "dark gas" needed to explain the 
relatively strong [CII] in diffuse clouds, we extended the ex- 



citation model of N(CII) versus I(CII) (equation 2) to a two- 
layer cloud, containing HI and H2. In the optically thin regime, 
which as shown above is a reasonable assumption for the dif- 
fuse sources, the total [CII] intensity is the sum of that in each 
layer, I(CII)totd= I(CII)hi + I(CII)h,. We can solve for I(CI)hi as 
follows: use A^(HI) = 1.82xlO''*I(HI)(cm-2) to solve for N(HI) 
from I(HI), then N(C+) = X(C+) N(HI); and, finally I(CII)hi can 
be determined from equation (2) assuming a temperature and 
density, (T<:,„, n(HI)) in the HI layer Since, I(CII)h, = I(CII)totai 

- I(CII)hi, we can now solve for N(C^)h,, using equation (2), as 
a function of temperature and density (TA^,„,n(H2)). For illustra- 
tive purposes in Figure 2, we have chosen to consider N(H2) a 
variable, and define the relative column densities of molecular 
and atomic gas through a parameter, a - N(H2)/N(HI). 

In Figure 2 we plot I(CII) versus I(HI) for different values 
of a (dashed lines) ranging from to 6 (a = corresponds to 
a pure HI cloud), for a fixed temperature and density (T,n(HI)) 

- (lOOK, 300 cm"-*), and assume constant pressure and temper- 
ature, so n(H2) = n(HI). We assume that I(HI) yields the col- 
umn density N(C^) in the HI region, and N(H2) is proportional 
to N(HI). A modest layer of H2 (N(H2) = 1 to 2 x N(HI)) fits 
the distribution of many of the stronger [CII] clouds. However, 
a much thicker layer of H2 (a = 4 to 6) is required to explain the 
strongest [CII] emitting clouds in this sample. 

To characterize and compare the clouds and their dark 
gas content, we calculated the visual extinction Ay and ra- 
tio N(H2)/N(HI) for each source, assuming a fixed (T, n) = 
(lOOK, 300 cm"^); these are plotted in Figure 3 for each source. 
We follow the procedure described above using N(HI) oc I(HI) 
to derive N(C^)hi, then calculate the corresponding I(CII)hi; 
subtract from I(CII)totai; finally, derive N(C^)h, using equa- 
tion 2, and N(H2) = N(C+)h,/(2X(C+)). then Ay = (N(HI) 
+2N(H2))x5.35xl0"2l The extinction ranges from 0.15 to 1.3 
mag. and about 75% of the sources have N(H2)> 0.5N(HI), and 
about 30% have N(H2)>2N(HI). If the density and/or tempera- 
ture in the H2 layer is lower than assumed here, N(H2) will be 
higher; and for higher (n,T) it will be lower. Here we assumed 
typical diffuse cloud conditions. However, [CII] preferentially 
identifies warm gas, so some sources may be hotter and denser, 
perhaps in proximity to UV sources or sampling shocked gas, 
and have smaller N(H2) than we calculate. Although it may the 
the case in a few individual clouds, statistically our sample of 
diffuse clouds observed in [CII] emission provides evidence for 
warm molecular "dark gas" in diffuse regions. 

4. Discussion 

The suggestion of the presence of "dark gas" in the interstellar 
medium is not new, Bloemen et al. (1986) used COS-B, Strong 
et al. (1996) EGRET, and Abdo et al. (2010) FERMI-LAT y-ray 
data to show that the Galaxy has more gas mass than indicated 
by HI and CO surveys alone. Joncas et al. (1992) compared HI 
and 100 /im dust emission in infrared cirrus and from the excess 
infrared emission suggested there might be molecular hydrogen 
or cold atomic gas (T < 30K) - see also Reach et al. (1994), who 
include CO. Grenier et al. (2005), who applied the term "dark 
gas" to this component, extended this approach to include y-ray 
observations and concluded that there were "dark gas" layers (H2 
without CO) surrounding all the nearby CO clouds. In addition, 
models of diffuse clouds and cloud envelopes (PDRs) predict 
significant amounts of H2 gas with no CO, but containing C^ 
and some neutral carbon, C° . 

Here we detected the warm "dark gas" with [CII], and 
showed it is a tracer of this important ISM component. Velusamy 
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Fig. 2. Observed [CII] intensities (filled circles) versus I(HI) 
(lower axis); the corresponding N(HI) is given on the upper axis 
for optically thin lines. Model calculations of I(CII) versus I(HI) 
as a function of (Ti.,„,n(HI),Q'), ■where a is the ratio of N(H2) to 
N(HI) in the cloud (see text). The labels (T,n,a) for each curve 
are given in the box. The solid lines are for clouds in which all 
hydrogen is in HI; the bottom line represents typical diffuse HI 
conditions, while the top solid line is an upper limit for atomic 
clouds. The dashed lines include H2 gas in addition to HI, (see 
text); with values of a'=N(H2)/N(HI) from 1 to 6. Only about 
50% of the sources fit within the limits defined by the purely HI 
models, the remaining sources require significant N(H2) to ex- 
plain the large [CII] emission. Most can be explained by N(H2) 
equal to (1 - 4)xN(HI). 



et al. (2010) use the GOT C+ data base to analyze [CII] asso- 
ciated with '^CO clouds and find evidence for warm H2 layers 
sandwiched between the HI and CO regions. Thus, warm H2 gas 
is an important component in interstellar clouds, in agreement 
with cloud-chemical models. We estimate the [CII] emission 
from all diffuse components (this paper and transition clouds 
Velusamy et al. 2010 and find their sum to be about equal to 
that of the dense PDRs Pineda et al. 2010) along these 16 LOS. 

5. Conclusions 

We have detected, for the first time, warm "dark gas" using the 
C^ ^P3/2-^Pi/2 fine structure emission line at 1.9 THz (158 yum). 
In just 16 lines-of-sight we detected 146 cloud components in 
[CII] emission. Of these, 29 are diffuse clouds, and more than 
75% likely contain warm H2. We analyzed the H2 content and 
found that it ranges from a fraction of the HI cloud component, 
to being the dominant component, with N(H2) = (1 to 4)xN(HI). 
If Tt,„ and n(H2) are less than assumed in Figure 3, N(H2) will be 
a larger fraction of the hydrogen in the cloud, and conversely for 
higher (n,T). Heterodyne observations of [CII] are an excellent 
probe for studying the warm "dark gas"; they have the advan- 
tage of locating the gas in the Galaxy, and providing excitation 
and dynamical information about the clouds. High spectral res- 
olution observations will help us understand the relative contri- 
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bution of [CII] from the diffuse ISM compared to dense PDRs 
associated with star forming regions. However, [CII] emission 
requires T^,,, > 30K, and will miss cold HI and H2 clouds. In 
conclusion, the warm "dark gas", is no longer "dark", indeed it 
shines quite brightly in C^. The full implementation of GOT C+ 
will detect thousands of such clouds, and yield information for 
understanding how their properties depend on ISM conditions. 
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